AD/ A- 000  06  4 


PRELIMINARY  REPORT  ON  AND  PROPOSAL  FOR 
FUTURE  FOCAL  MECHANISM  STUDY  OF  NORTH 
ATLANTIC  EARTHQUAKES  BY  THE  SURFACE 
WAVE  METHOD 

K . A k i , ct  a 1 

Royal  Norwegian  Council  for  Scientific  and 
Industrial  Research 


P r e p a r e d for: 

Massachusetts  Institute  of  Technology 
Air  Force  Technical  Applications  Center 
Advanced  Research  Projects  Agency 


10  September  1974 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


REPORT  DOCUMENTATION  PAGE 


report  number 


2 GOVT  ACCESSION  NO. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3 RECIPIENT’S  CATALOG  NUMBER 


FO  86  06-74-000  49 


4 TITLE  (and  Subtitle 


s type  of  report  a period  covered 


Preliminary  report  on  and  proposal 
for  future  focal  mechanism  study  of  North 
Atlantic  earthquakes  by  the  surface  wave 
me thod . 


One-time  technical 


6 PERFORMING  ORG  REPORT  NUMBER 

Scientific  Rep.  2-74/75 


A j T H O R i 


0 CONTRACT  or  grant  NUMBERfs; 


K.  Aki,  M.I.T.,  Cambridge,  Mass.,  U.S.A. 
E.S.  Husebye,  NTNF/NORSAR,  Kjeller,  Norw, 


F08606-74-C-0049 


9.  PERFORM  N G ORGANIZATION  NAME  AND  ADDRESS 


NTNF/NORSAR 
Post  Box  51 
Kjeller,  Norway 


11  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

VELA  Seismological  Center 
312  Montgomery  Street 

Alexandria,  Va.  22314  


14  MONITORING  AGENCY  NAME  0 ADDRESS  //  Jifterent  from  Controlling  Office) 


10  PROGRAM  ELEMENT  PROJECT,  TASK 
AREA  0 WORK  UNIT  NUMBERS 


NORSAP  Phase  3 


I?  REPORT  DATE 


10  September  1974 


15.  SECURITY  Class  (of  this  report) 


1 5a,  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


17  DISTRIBUTION  STATEMENT  (of  the  mhstrmct  entered  in  Block  20,  It  different  from  Report) 


19  KEY  WORDS  'Continue  on  reverse  side  if  nece  ssary  mnd  Identify  by  block  number) 


Seismicity  and  tectonics  of  the  North  Atlantic,  Focal 
mechanisms  and  depths  of  small  earthquakes  from  surface  wave. 
Rayleigh  and  Love  wave  dispersion  and  attenuation. 


20  ABSTRACT  f Continue  on  reverse  side  If  necessary  and  identify  hv  block  number) 


We  are  interested 
the  Atlantic  Ocean 
in  number  of  analy 
gion  where  >the  sei 
it  is  demonstrated 
small  earthquakes 
as  reference  event 


in  the  seismicity  and  earthquake  mechanism  of 
north  of  Iceland.  In  this  respect  an  increase 
zable  earthquakes  is  badly  needed  for  this  re- 
smic  activity  is  relatively  low.  In  this  report 
how  the  surface  wave  method  can  be  extended  to 
using  known  fault  plane  solutions  in  the  region 
s.  The  outcome  of  these  studies  - Rayleigh  and 


DD  ,:°NBMT3  1473  EDITION  OF  1 NOV  6S  IS  OBSOLETE  ' 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fHTien  Data  Entered) 


NATION  A!  TECHNIC  A! 
INFORMA I ION  SERVIff 


Siv  u i/bt  I ' VA 


J 


UCUWITY  CLASSIFICATION  OF  Thu  PAQt(Whmt  Df  Ani.r.rfj 


ii 


AFTAC  Project  Authorization  No.: 

ARPA  Order  No.  : 

Program  Code  No.  : 

Name  of  Contractor  : 

Effective  Date  of  Contract 
Contract  Expiration  Date 
Contract  No. 

Project  Manager 
Title  of  Work 

Amount  of  Contract  : 

Contract  period  covered  by 

the  report  : 


VT/5702/B/ETR 

2551 

5F10 

Royal  Norwegian  Council 

for  Scientific  and  Industrial 

Research 

1 July  1974 

30  June  1975 

F08606-74-C-0049 

Nils  Mar^s  (02)  71  69  15 

Norwegian  Seismic  Array 
(NORSAR)  Phase  3 

$900  000 

1 July  1974-30  June  1975 


The  views  and  conclusions  contained  in  this  document  are  those 
of  the  authors  and  should  not  be  interpreted  as  necessarily 
representing  the  official  policies,  either  expressed  or  im- 
plied, of  the  Advanced  Research  Projects  Agency,  the  Air 
Force  Technical  Applications  Center,  or  the  U.S.  Government. 


This  research  was  supported  by  the  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense  and  was  monitored  by 
AFTAC/VSC,  Patrick  AFB  FL  32925,  under  Contract  No.  F08606- 
74-C-0049 . 


iii 


CONTENTS 

Page 

PROBLEMS  AND  HYPOTHESES  1 

METHOD  5 

USE  OF  SURFACE  WAVES  5 

DATA  7 

PRELIMINARY  RESULTS  7 

(1)  The  Mohn's  ridge  earthquake  of 

31  May  1971  - Event  1 7 

(2)  Amplitude  spectrum  of  Rayleigh  v/aves  10 

(3)  Phase  spectrum  of  Rayleigh  waves  13 

(4)  Amplitude  spectrum  of  Love  waves  17 

(5)  Phase  spectrum  of  Love  waves  20 

(6)  Instrument  response  for  VLP 

seismograph  at  KON  22 

FUTURE  STUDIES  24 

REMARK  25 

ACKNOWLEDGEMENT  26 

REFERENCES  27 


PROBLEMS  AMD  HYPOTHESES 


We  are  interested  m the  seismicity  and  earthquake  mechanism 
of  the  Atlantic  Ocean  north  of  Iceland.  Fig.  1 shows  the 
nap  of  epicenters  determined  for  528  earthquakes  during  the 
period  1955-1972,  reproduced  from  Husebye  ot  al  (in  press). 
They  are  shown  in  relation  to  major  tectonic  features  such 
as  the  mid-ocean  ridges,  fracture  zones  and  shel^  edqes. 

We  are  particularly  interested  in  the  mechanism  of  intra- 
plate earthquakes  occurring  in  the  Norwegian  Sea  between 
the  Mohn's  ridge  and  the  Norwegian  const. 

The  mechanism  of  these  intra-plate  earthquakes  was  discussed 
by  Husebye  et  al  in  relation  to  geological  and  geophysical 
data  from  the  area  within  the  context  of  plate  tectonics. 

Their  discussions  suggest  the  following  alternative  hypotheses 
for  their  mechanism. 

1)  The  post-opening  marginal  subsidence;  in  this  case, 

we  expect  a normal  faulting  with  fa-lt  strikes  parallel 
to  the  coast. 

2)  The  intra-plate  stress  which  appears  to  have  the  hori- 
zontal compressive  axis  parallel  to  the  direction  of 
plate  motion,  as  pointed  out  first  by  Mendiguren  (1971) 
for  the  Nazca  plate  and  later  elaborated  by  Sykes 

and  Sbar  (1973).  In  this  case,  we  expect  a thrust 
faulting  with  fault  strikes  perpendicular  to  the 
direction  of  plate  motion  or  a strike-slip  faulting 
with  the  pressure  axis  parallel  to  the  direction  of 
plate  motion. 
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Fiq.  1 Epicenter  distribution  (1955-1972)  in  the  Norwegian-Greenland 
Sea  and  adjacent  areas.  Epicenters  and  selected  sea  bottom 
contours  (in  meters)  are  shown  together  with  main  structural 
features.  Earthquakes  for  which  focal  mechanisms  have  been 
published  are  numbered  from  1 to  15  (Husebye  et  al,  in  press). 
The  two  earthquakes  used  in  analysis  in  this  study  are  marked 
by  crosses. 
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3)  The  seismic  zone  in  the  Norwegian  Sea,  especially  the 
belt  v/hich  continues  southward  as  an  apparent  exten- 
sion of  the  Knipovich  ridge,  may  be  a plate  houndary. 
As  shown  schematically  in  Fig.  2,  if  the  mid-Arctic 
ridge  spreads  faster  than  the  Mohn's  ridge,  we  expect 
right-lateral  strike  slip  along  the  Knipovich  ridne 
as  well  as  along  its  extension  into  the  Norwegian 
Sea.  Two  fault  plane  solutions  published  by  Lazareva 
et  al  (1965)  in  this  zone  are  right-lateral  strike- 
slip  along  fault  with  the  North-South  strike,  con- 
sistent with  this  hypothesis. 

4)  If  the  Mohn’s  ridge  is  spreading  faster  than  the  mid- 
Arctic  ridge,  v/e  expect  left-lateral  strike-slip 
along  the  extension  of  the  Knipovich  ridge. 

^n  order  to  examine  these  hypotheses,  we  shall  determine 
the  focal  mechanism,  focal  depth  and  other  source  para- 
meters as  well  as  the  surface  wave  dispersion  and  attenua- 
tion properties  in  various  parts  of  the  North  Atlantic. 


METHOD 


As  well  known,  reliable  fault  plane  solutions  can  be  ob- 
tained only  for  relatively  large  earthquakes  (magnitude  > 

5.5)  which  give  clear  onset  of  P waves  on  WWSS  LP  records. 

Our  first  step  is  to  find  as  many  a-  possible  earthquakes 
with  ’-r.cwn  fault  plane  solutions  in  the  region  of  our 
interest.  We  shall  use  them  as  reference  events  when 
applying  the  surface  wave  method  to  smaller  earthouakes. 

The  surface  wave  method  can  give  us  an  accurate  determination 
of  focal  depth  of  an  earthquake  if  the  fault  plane  solution 
is  known  (Weidner  and  Aki,  1973).  We  can  also  determine 
the  phase  velocity  and  attenuation  of  surface  waves  along 
the  path  from  an  epicenter  to  a station  if  the  fault  plane 
solution  is  known.  Once  the  phase  velocity  and  attenuation 
are  known  for  various  parts  in  the  region  of  our  interest, 
we  can  determine  the  source  factor  of  amplitude  and  phase 
spectra  for  any  earthquake  in  the  area.  If  we  have  both 
Love  and  Rayleigh  wave  source  spectra  for  a few  different 
azimuths,  we  can  probably  make  a reliable  determination 
of  the  focal  mechanism,  seismic  moment  and  focal  depth. 
Additional  data  from  body  waves,  such  as  the  pP-P  time 
difference  and  the  sense  of  first  P motion  at  crucial 
directions,  would  be  helpful  to  increase  the  accuracy 
of  the  determinations.  The  outcome  of  the  above  studies  - 
the  dispersion,  attenuation,  focal  depth,  seismic  moment 
and  focal  mechanism  - are  useful  quantities  for  discussing 
the  tectonics  of  the  North  Atlantic. 


USE  OF  SURFACE  WAVES 

In  analyzing  surface  waves,  we  shall  assume  that  an  earth- 
quake is  a slip-dislocation,  and  the  earth  is  a laterally 
homogeneous  elastic  medium.  For  periods  longer  than  10 


<5 


seconds  and  for  earthquakes  with  M < 6 , we  shall  further 

assume  that  the  fault  slip  is  a step-function  an  time  and 

effect  of  finite  size  of  fault  nav  be  noqloctod  (Tsai 

and  Aki,  1970).  We  shall  use  a layered  medium  aonropri ate 

for  the  North  Atlantic.  The  model  used  in  calculations  so 

far,  denoted  the  Norwegian  Deep  Sea  Basin  model,  is  listed 
in  Table  1 


NORWEGIAN 

DEEP  SFA 

basin  model 

Layer 

Thickness 

(km) 

Density 

(g/cm3) 

P-Veloci  ty 
(km/sec) 

S-Veloci ty 
(km/  sec) 

1 

3.50 

1.00 

1.50 

0.0 

2 

1.30 

2.10 

2.10 

1.05 

3 

1.50 

2.40 

4.60 

2.40 

•1 

4.00 

3.09 

7.40 

4.28 

5 

26.7 

3.40 

8.  11 

4.61 

6 

20.0 

3.  37 

8.01 

4.56 

7 

20.0 

3.  37 

7.95 

4.56 

8 

29.0 

3.  37 

7.71 

4.40 

9 

20.0 

3. 30 

7.68 

4.  34 

TA3LF  1 


Crust  and  upper  mantle  model  for  the  Norwegian  Sea  used  in 
■is  study.  For  Love  wave  calculations,  the  water  layer  is 

^nzVU972rindCF^n^ry  CrUStal  layerS  are  modeled  after 

density  Li  s w , ‘ ?nd  Windish  {in  . Note  that  the 

ensiti  and  S-wave  velocity  of  the  sedimons  represent  an 

lligent  guess.  The  upper  mantle  is  tvpical  for  oceanic 
regions  (Anderson,  1964).  oceanic 


Two  computer  programs  have  been  written  at  Massachusetts 
Institute  of  Technology  (M.I.T.)  for  calculating  the  ampl 
tude  and  phase  spectra  of  both  Love  and  Rayleigh  waves  fo 
the  above  models  of  source  and  medium.  The  first  one  is 
due  to  M.  Salto,  and  computes  the  phase  velocities  and 
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eigenfunctions  for  an  arbitrary  layered  medium.  The  second 
one  is  written  by  Y.D.  Tsai,  and  computes  the  amplitude 
and  phase  spectra  of  Love  and  Ravleigh  waves  for  an  arbi- 
trary dislocation  source  using  the  output  of  the  Saito 
program.  Tsai's  program  is  based  on  the  theory  of  Saito 
(1967) , which  is  explained  in  detail  in  a class  note  by  K. 
Aki.  Both  programs  are  working  satisfactorily  on  the 
computer  at  NORSAR. 


DATA 

Husebye  et  al  (in  press)  were  able  to  collect  only  15  fault 
plane  solutions  from  about  15  years'  earthquake  data  for 
the  area.  We  are  trying  to  obtain  additional  fault  plane 
solutions  by  collecting  and  reading  the  'fWSS  and  other 
records  by  ourselves.  Reliable  fault  plane  solutions  are, 
however,  expected  only  for  relatively  large  magnitudes 
(mb  ~ 5.5).  On  the  other  hand,  the  WNSS  LP  seismographs 

in  Iceland,  Greenland,  Spitsbergen  and  Fennoscandia  can 
register  long-period  Rayleigh  and  Love  waves  fron  earth- 
quakes with  Mg  around  5.  Further,  the  VLP  seismograph 
at  KON  can  record  long  period  waves  from  earthquakes  with 
Ms  around  4.  Thus,  the  surface  wave  method  can  be  ex- 
tended to  earthquakes  with  magnitude  one  or  two  units  lower 
than  the  P— wave  fault  plane  method,  thereby  increasing  the 
number  of  analyzable  earthquakes  by  more  than  an  order  of 
magnitude.  The  increase  in  number  of  analvzable  earthquakes 
is  badly  needed  for  the  North  Atlantic  region,  where  the 
seismic  activity  is  relatively  low. 


PRELIMINARY  RESULTS 

(1)  The  Mohn's  Ridge  Earthquake  of  31  May  1971  - Event  1 

According  to  Conant  (1973) , the  fault  plane  solution  of 
this  earthquake  consists  of  two  nodal  planes ) one  striking 
N42  E with  dip  angle  54°,  the  other  striking  N52°E  with 
dip  angle  64°.  The.  focal  parameters  for  this  event  are  given 
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in  Table  2.  If  the  former  is  the  fault  plane,  the  fault 
is  a normal  fault  with  small  right-lateral  strike-slip. 
The  Tsai  program  requires  the  slip  angle  b,  which  is 


EVENT 

DATE 

(m/  d/  y) 

0.  TIME 
(h  /m  /sec) 

LAT 

(deg) 

LONG 

(deg) 

DEPTH 

(km) 

MAG 

(n^) (Ms) 

REGION 

1 

05/31/71 

03. 4€ .52.2 

72. 2N 

1.2E 

33 

5.5  5.7 

Mohn • s 
Fiqe 

2 

05/21/72 

07.57.05.2 

73. 6N 

7.6F 

33 

4.5  - 

Mohn- 
Knipovich 
Ri  dge 

TABLE  2 

NOAA  focal  parameters  for  the  two  earthquakes  used  in  analysis 


related  to  the  difference  AY  between  the  strikes  of  two 
nodal  planes  and  the  dip  angle  6 of  the  auxiliary  plane 
(when  one  of  the  nodal  planes  is  called  "fault  pl«ne" , the 
other  is  called  "auxiliary"}  by  the  following  formula: 

cos  <(>  = sin  AY  • sin  5 

For  Event  1,  AY  = ]0°,  6 r.  64°,  theref  -e  cos  <j> 

* = 81°. 

The  output  of  Tsai's  program  gives  the  amplitude  and  phase 
spectra  at  a certain  azimuth  measured  counterclockwise 
from  the  strike  of  fault  plane. 

We  choose  the  strike  direction  X as  the  positive  X-axis, 
the  upward  vertical  as  the  positive  Z-axis.  The  positive 
Y-axis  in  the  right-hand  system,  then,  is  fixed  in  the 
space.  We  define  the  dip  angle  as  the  angle  between  the 
fault-dip  direction  and  the  positive  Y-axis.  For  the  Mohn's 
ridge  earthquake,  taking  the  X-axis  to  N42°E , the  positive 
Y-axis  is  toward  N48°W.  Since  the  fault  plane  is  dipping 
to  NW,  the  dip  angle  is  54°. 


(1) 

= 0.156, 
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Fig.  3 Comparison  of  observed  and  theoretical  Rayleigh  wave  spectra. 

Thi.  Rayleigh  waves  were  recorded  at  KON  for  the  Mohn's  ridge 
earthquake  of  31  May  1971.  The  fault  plane  solution  of 
Conant  (1973)  was  used  .n  the  theoretical  calculations, 
assuming  a focal  depth  of  8 km.  Fault  plane  strike  is 
N42  E with  dip  angle  54  , auxiliary  plane  strike  is  N52  E 
with  dip  angle  64  . 


The  slip  angle  is  the  angle  between  the  X-axis  and  the 
direction  of  slip  of  the  body  containing  the  positive  Y- 
axis  relative  to  the  body  containing  the  negative  Y-axis. 

The  angle  is  measured  from  the  positive  X-axis  counter- 
clockwise in  the  fault  plane.  The  slip  angle  for  the  Mohn's 
ridge  earthquake  is  279°. 

Tsai's  program  lists  the  amplitude  | A ( to)  | and  source  phase 
delay  <J)(u>)  under  the  heading  of  POINT  AMP.  and  FOCAL  PHASE 
respectively.  For  example,  for  Kongsberg  A=*2000  km  and 
AZIMUTH=*240°  from  the  Mohn's  ridge  earthquake  which  is 
taken  at  a focal  depth  of  5 km  (below  sea  level) , 
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A(u)  0.835  x 10  26  cmsec/dynecm  and  <Mw)  =0.377  at  frequency 
0.02  Hz.  A(u)|  is  the  di splaconont  amplitude  spectral 
density  corresponding  to  a unit  moment  earthquake,  *(U) 
is  measuie  in  parts  of  a cycle.  They  are  defined  as 

A(t)  TT  f |A(ul)  oxp  ( iut  - i - i ( 0} ) ) dui  (2) 

where  A(t)  is  the  ground  displacement  in  cm  (already 

corrected  for  instrument  response)  and  C((„)  is  the  phase 
velocity . 


■iii Amplitude  Spectrum  of  Rayleigh  Waves 

Fig.  4 shows  the  amplitude  spectral  density  |A(w)|  for 
Rayleigh  waves  from  the  Mohn’s  ridge  earthquake  recorded 


Fig.  1 


Dependence  of  amplitude  spectrum  of  Rayleigh  waves  on 
focal  depth.  The  model  is  the  Mohn’s  ridge  earthquake 
recorded  at  Kongsberg  (KON) . S 


11 


at  KON.  The  seismogram  trace  f(t)  is  digitized  and 
A(w)  is  obtained  as 

A(u)  = I / f(t)elu:t  dt  | / I(u)|  (3) 

where  I(w)  is  the  magnification  of  the  seisnonraph  at 
frequency  to,  which  is  the  ratio  of  record  amplitude  to 
ground  displacement. 

The  observed  spectrum  is  compared  in  Fig.  3 with  the 
theoretical  one  calculated  for  the  fault  Diane  solution 
given  above  using  the  Tsai  program.  The  theoretical  curve 
is  given  for  a source  with  the  seismic  moment  1 x 1025  dyne  cm 
buried  at  a depth  (below  sea  level)  of  0 km. 

The  estimation  of  seismic  moment  depends  on  the  assumption 
on  (1)  focal  depth,  (2)  focal  mechanism,  especially  the 
azimuth  of  fault  stride , and  (3)  attenuation.  The  effect  (3) 
may  be  neglected  for  periods  longer  than  20  sec  for  the 
epicentral  distances  we  are  interested  in  here.  The  effect 
(1)  is  very  complicated  as  shov/n  in  Fig.  4.  This  effect 
requires  a very  careful  examination.  Fig.  4 suggests, 
for  example,  the  best  fit  between  theoretical  and  observed 
may  be  obtained  for  depths  greater  than  55  km.  This 
was,  as  a matter  of  fact,  a conclusion  given  in  Tsai's 
thesis  (1969)  that  the  dip-slip  earthquakes  in  the  mid- 
ocean ridges  are  deep.  Later,  Weidner  and  Aki  (1973) 
showed  by  the  combined  use  of  phase  and  amplitude  spectra 
that  these  earthquakes  must  he  shallow.  The  discrepancy 
at  T<20  sec  between  the  observed  and  theoretical  amplitude 
such  as  shov/n  in  Fig.  4 is  attributed  to  the  attenuation 
due  to  soft  sediment  by  Weidner  (1972). 
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Fig.  5 Rayleigh  wave  radiation  from  the  Mohn’s  ridge 

earthquake.  For  details  on  focal  parameters  see 
the  caption  of  Fig.  3. 


The  effect  (2)  can  be  examined  with  the  use  of  Fig.  5, 
which  shows  the  radiation  pattern  of  |a(u>)  | for  all  azimuths. 
A small  change  in  the  strike  of  the  fault  which  is  not  very 
well  constrained  for  this  earthquake  can  cause  a large 
change  in  |A(w) | near  the  nodal  direction.  To  resolve 
the  above  ambiguities  due  to  uncertain  focal  mechanism 
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the  phase  spectra  of  Rayleigh  waves  as  well  as  the  ampli- 
tude and  phase  spectra  cf  hove  waves. 


(3)  Phase  Spectrum  of  Rayleigh  Waves 

The  phase  spectrum  of  Rayleigh  waves  is  a very  nowerful 
tool  for  determining  the  focal  depth.  For  examole,  in  the 
case  of  Mohn's  ridge  earthquake,  the  theoretical  source 
phase  at  50  seconds  period  is  roughlv  0.375  for  all  azi- 
muth for  depths  h<25  km,  and  changes  to  0.7  ~ 0.8  for 
depths  h>35  km.  This  difference  can  be  easily  detected 
once  reliaole  dispersion  curve  for  Rayleigh  waves  is 
established  for  the  area. 

From  the  record  f(t)  of  Rayleigh  waves  (vertical  component, 
upward  positive)  we  find  the  phase  delay  £ (...)  , defined  as 

f(t)  = / !f(u)  eiu,t"ivSU’)  du  (4) 

where  t is  measured  from  the  origin  time  of  an  earthquake. 

For  a single  mode,  we  can  write 

*U)  = - * V-')  + ' 2N'  <5) 

where  x is  the  epicentral  distance,  C(-)  is  the  phase  velocity, 
: INST  ^ ““  ^ *s  the  PhasG  delay  of  record  trace  relative  to 
ground  displacement.  For  the  WWSS  LP  record, 
approaches  v/2  as  (This  can  be  used  to  check  the 

" error  in  phase  calculations.) 

In  the  output  of  Tsai's  program,  £c(^)  is  called  the  FOCAL 
PHASE.  For  simple  fault  geometries,  the  values  of  itQ ( u>) 
are  shown  in  Figs.  6 and  7.  They  can  be  used  to  check 
blunders  such  as  the  tt  error. 
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SOL  RCE  PHASES  FOR  SIMPLE  FAULT  GEO'TETRICS 

a l_Jupward_postive]__  Ray  lei  ghjyave_  displacement 
AR(t)  = / lAR(u)  |exp[iu)t-iu)A/CR(u>)-$R]  du> 

Ir3DSy®£se__(counterclockwise_positiye)__Loye_waye_displacenient 

AL(t)  = J exp[iut-iuA/CT  (w)-a  } du» 

L 


(1)  Strike-slip  on  a vertical 


fault 


iq.  b Sourco  phases  for  Love  and  Rayleigh  waves  for  simple 

fault  geometries,  i.e.,  strike-slip  on  a vertical  fault 
u *S  GP1  central  distance,  C<ui)  and  4.  are  phase  velocity’ 
and  phase  delay  for  respectively  Rayleigh  or  Love  waves. 


Fron  the  observed  phase  * (u) , one  can  obtain  the  phase 
velocity  CU)  if  the  source  phase  :q(u)  is  known.  For 
Event  1,  we  assumed  that  the  focal  depth  is  less  than 
25  kn  and  therefore  $o(u>)  = 0.375  at  T=*50  sec.  For 
different  choices  of  integer  N,  one  gets  the  values  of 

Clu)  at  T=50  sec*'  such  as...,  3.51,  3.99,  4.62,  ...  . 
Obviously,  3.51  km/sec  is  too  low  and  4.62  is  too  high. 
They  are  both  definitely  unacceptable.  This  fixes  the 
integer  N,  and  determines  uniquely  the  entire  phase- 
velocity  curve  as  shown  in  Fig.  8. 


The  observed  phase  velocity  curve  is  compared  with  the 
theoretical  one  used  for  the  eigenfunction  calculation 
by  the  Saito  program.  Clearly,  the  path  from  the  Mohn’s 
ridge  to  KON  is  a mixed  one  (o.g.,  see  Talwani  and  Eldholm, 
1972)  and  the  observed  curve  departs  from  the  theoretical 
one  for  the  ocean  considerably  at  short  period. 


SOURCE  PHASES  FOR  SIMPLE  FAULT  GEOMETRICS 


(2)  Dip-slip  on  an  inclined  fault 


THRUST  t-AULl 


Ei  7*  Source  phases  for  Love  and  Rayleigh  waves  f r simple 

fault  geometries,  i.e.,  dip-slip  on  an  inclined  tault. 


Suppose  that  the  focal  depth  was  not  less  than  25  km 
but  greater  than  35  km.  Then,  the  source  phase  changes 
by  about  ",  giving  alternative  values  of  phase  velocities 
3.73  and  4.28  at  T~I 0 sec.  These  values  seen  to  be 
rather  too  low  or  too  high,  but  additional  data  will  be 
needed  to  resolve  this  definitely. 


MOWAt . AN  it  A OtlP  BAb<N 


CANADIAN  V«KD 


'SStHvED-evfM  ' 


RAU  t ->»  *A«E  PS^R^ON 


■ PERIOD  I sec) 


Rayleigh  wave  phase  velocity  obtained  by  analysis  of 
the  wi in ' s ridge  earthquake  of  31  May  1971  (Event  1) 
usinc,  Kongsberg  (KON)  records.  The  Norwegian  Sea  Deep 
Basin  dispersion  curve  is  calculated  from  the  model  in 
Table  1,  while  the  Canadian  Shield  curve  is  taken  from 
Brune  and  Dornan  (1963). 


AMPilTuCN  I AJ <»  ! i . ' <J »*>***1 
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(4)  Anplitudo  Spectrum  of  Love  Waves 

The  anplitudo  spectrum  of  Love  waves  does  not  depend 
critically  on  the  focal  depth  as  Rayleiqh  waves  do. 
Fig.  9 shows  the  dependence  for  the  case  of  record 

for  the  Mohn's  ridge  earthquake. 


Fig.  9 Love  wave  amplitude  dependence  on  focal  depth  for  the 
Mohn's  ridge  earthquake  recorded  at  Kongsberg. 


The  observed  spectrum  is  compared  with  the  theoretical  in 
Fig.  10.  A good  fit  of  the  amplitude  spectral  shape  is 
obtained  for  20  < T * 50  sec.  The  observed  spectral 
density  at  T * 50  sec  is  around  0.35  • 10"*1  cm  sec,  and 
the  theoretical  density  is  0.17  • 10  cn  sec/dyne  cm. 

This  gives  a seismic  moment  of  =*  2 • 10^  dyne  cm. 

This  is  about  a factor  of  2 grocit^r  i$an  the  value  previously 
estimated  from  Rayleigh  waves. 


AMF'LITIj[H  i A(uil  crr>  s**c 
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Fiq.  IQ  Observed  Love  wave  amplitude  spectrum  for  the  Mohn’s 
ridge  earthquake  of  31  May  1971,  using  the  Kongsberg 
sta  taon. 


The  above  discrepancy  nay  be  resolved  if  we  allow  a small 

change  in  the  fault  plane  solution.  As  shown  in  Fig.  11 

and  also  in  Fig.  5 both  Love  and  Rayleigh  wave  amplitudes 

are  sensitive  to  the  change  ir.  the  azimuth  of  fault  strike. 

For  example,  if  we  rotate  the  fault  clockwise  by  20°  or  so, 

making  the  strike  direction  a*-  N20°E,  then  both  Rayleigh 

and  Love  waves  would  give  the  same  moment  estimate,  about 
25 

1.5  • 10  dyne  cm.  (Since  the  distance  dependence  of 
surface  wave  amplitude  is  inversely  proportional  to  the 
square  root  of  distance,  and  the  theoretical  value  is 
calculated  for  A=2000  km,  but  the  actual  distance  is 
about  1500,  the  moment  is  overestimated  by  J ~~  = 1.15.) 
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(5)  Phase  Spectrum  of  Love  Waves 

For  Love  waves,  we  define  the  focal  phase  ^(w)  in  the 
same  way  as  for  Rayleigh  waves.  The  Lo^e  wave  displacement 
f(t)  (horizontal  transverse  component,  counterclockwise 
as  seen  from  the  source  positive)  is  given  by 

f(t)  = / F(u)  exp<  iut-iMw)  } dw  (6) 


where 


:(u)  = 


rrr 


(w) 


+ 'ISST1'1’  + 2“” 


(7) 


For  the  Mohn's  ridge  earthquake,  £0(oi)  is  +0.13  cycle  at 
T^50  sec.  Different  choices  of  integer  N give  values  of 
the  phase  velocity  C(w)  as  ...,3.83,  4.40,  6.20,...  . Only 
possible  value  is  4.40,  and  this  choice  determines  the 
rest  of  the  phase  velocity  curve,  as  shown  in  Fig.  12. 


The  phase  velocities  for  periods  less  than  20  sec  show 
some  irregular  behavior.  The  contamination  by  higher 
nodes  or  Rayleigh  waves  nay  be  effective  in  the  short 
periods . 

In  addition  to  the  Mohn's  ridge  earthquake,  we  analyzed 
the  KON  record  of  a shock  at  the  intersection  of  Knipovich 
ridge  and  Mohn's  ridge  (Event  2)  which  occurred  on  21  May 
1972.  The  focal  parameters  of  this  event  are  given  in 
Table  2.  I f we  assume  that  the  source  mechanism  of  this 
earthquake  is  the  same  as  the  two  nearby  earthquakes  given 
by  Lazareva  et  al  (1965),  then  we  expect  that  the  focal 
phase  <t)Q  of  Love  waves  to  be  +.125  at  KON. 


PHASE  VELOCITY  ( km /sec  I 


nonwe  -'A*  se  a ot1  p w.in 


observed-eve.’  * 


OBSERVE  0 
Event  ; 


* •^utenuewl.  birch  n 
/ CONT  MENTAL  model 


LOVE  WAVE  DISPERSION 


PERIOD  (i»cl 


Fig.  12  Love  wave  phase  velocity  measurements  using  Kongsberg 
records  for  the  Mohn's  ridae  earthquake  of  31  May  1971 
(event  1)  and  Mohn-Knipovich  ridge  earthquake  of  21  May 
1972  (event  2) . The  Norwegian  Sea  Deep  Basin  dispersion 
curve  is  calculated  from  the  model  in  Table  1,  while 
the  Gutenberg-Birch  II  curve  is  taken  from  Anderson  (1964) 
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The  observed  phase  spectrum  of  Love  waves  from  Event  2 gives 
the  phase  velocity  of  4.47  km/sec  at  T=50  sec  for  <£s  = +.125. 
The  corresponding  phase  velocity  curve  is  also  shown  in 
Fig.  12. 

If  v:e  choose  the  hypothesis  (4),  and  assume  the  strike- 
slip  motion  opposite  to  those  given  by  Lazareva  et  al 
(1965)  (and  therefore  <}>  = +.625),  we  get  the  phase  velocity 

values  at  T«50  sec  such  as  4.10  and/or  4.82  km/sec.  They 
seem  to  be  too  low  or  too  high,  not  in  favor  of  the  hypo- 
thesis (4)  . 


(6)  Instrument  Response  for  VLP  Seismograph  at  KON 

Unfortunately,  the  phase  delay  response  of  the  VLP  seis- 
mograph at  KON  (Locke,  1971)  is  not  available  at  this 
moment.  We  tried  to  obtain  it  empirically  by  comparing 
the  records  of  VLP  seismographs  with  those  of  LP  seismo- 
graphs at  the  same  station.  We  digitized  and  Fourier- 
analyzed  both  Love  and  Rayleigh  waves  from  the  Mohn- 
Knipovich  ridge  earthquake,  and  measured  the  amplitude 
spectral  ratios  and  phase  differences  at  various  frequencies. 
Combining  these  measurements  with  the  calculated  response 

for  the  standard  WWSS  LP  (T  =15  sec,  e =1,  T =100  sec, 

o o g 

eg=0.93  and  coupling  coefficient  o=0.047  (Anonymous,  1966)), 
we  obtained  a reasonably  stable  estimate  of  amplification 
(Fig.  13)  and  phase  delay  curve  (Fig.  14)  for  the  VLP 
seismograph.  These  curves,  however,  should  be  considered 
as  tentative. 


SEISMOGRAPH  AMPLIFICATION 
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Fig.  14  Estimated  phase  delay  curves  for  the  Very  Long  Period 
(VLP)  seismograph  at  Kongsberg.  For  comparison  a 
similar  curve  for  the  ordinary  long  period  (LP)  .tation 
is  also  given. 

FUTURE  STUDIES 

The  above  preliminary  results  offer  the  first  approximations 
to  the  basic  quantities  involved  in  the  analysis  of  surface 
waves  from  the  North  Atlantic  ocean  for  the  purpose  of  source 
mechanism  studies.  We  made  a thorough  analysis  of  amplitude 
and  phase  spectra  of  Love  and  Rayleigh  waves  for  the  KON 
record  of  a Mohn's  ridge  earthquake  with  known  fault  plane 
solution.  From  this  analysis,  we  obtained  the  first  aoproxi- 
mation  to  the  phase  velocity  of  Love  and  Rayleigh  v/aves 
over  the  Norwegian  Sea. 


- 75  - 

The  next  step  will  be  the  following: 

1)  A similar  analysis  of  Rayleigh  and  Love  waves  from 
the  same  Mohn's  ridge  earthquake  recorded  at  other 
stations . 

2)  A similar  analysis  of  other  North  Atlantic  earthquakes 
with  known  fault  plane  solution  published  by  others 

or  obtained  by  ourselves.  Improve  the  phase  velocity 
curves  for  various  paths. 

3)  Determination  of  source  mechanism  of  smaller  earth- 
quakes from  the  measurement  of  source  phase  using  the 
first  approximation  phase  velocity  curves.  Combine 
Love  and  Rayleigh  waves  for  more  reliable  estimation. 

4)  Determination  of  seismic  moment  for  smaller  earthquakes 
using  the  source  mechanism  obtained  from  the  source 
phase  analysis.  Combine  Love  and  Rayleigh  waves  for 
more  reliable  estimation. 

REMARK 

The  value  of  seismic  moment  obtained  for  Event  1 earthquake 
is  comparable  to  the  Parkfield,  California,  earthquake 
for  which  detailed  source  studies  have  been  made  by  Aki  (1972). 
We  are  somewhat  bothered  by  this  result,  because  Mg  of 
Parkfield  earthquake  w as  determined  to  be  6.4  by 
Wu  (1968)  but  Mg  of  this  earthquake  seems  about  1 unit 
smaller  than  the  Parkfield.  We  could  not  find,  however, 
any  blunder  in  the  estimate  of  seismic  moment. 
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